Theory of cavity-polariton self-trapping and optical strain in polymer chains 
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We consider a semiconductor polymer chain coupled to a single electromagnetic mode in a cavity. 
The excitations of the chain have a mixed exciton-photon character and are described as polaritons. 
Polaritons are coupled to the lattice by the deformation potential interaction and can propagate 
in the chain. We find that the presence of optical excitation in the polymer induces strain on the 
lattice. We use a BCS variational wavefunction to calculate the chemical potential of the polaritons 
as a function of their density. We analyze first the case of a short chain with only two unit cells in 
order to check the validity of our variational approach. In the case of a long chain and for a strong 
coupling with the lattice, the system undergoes a phase transition corresponding to the self-trapping 
of polaritons. The role of the exciton spontaneous emission and cavity damping are discussed in the 
case of homogeneous optical lattice strain. 
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I. INTRODUCTION 

Recent advances in nano-spectroscopy have demon- 
strated the possibility of addressing a single conjugated 
polymer chain in a solid matrixii*^ These optical studies 
of single chains overcome limitations due to ensemble av- 
eraging and inhomogeneous broadening and give a clear 
picture of the exciton dynamics. Most interestingly, a ID 
singularity in the optical density of states, a T 1 / 2 tem- 
perature dependence of the the exciton lifetime^ and a 
macroscopic quantum spatial coherence^ have been ob- 
served in isolated red polydiacetylene chains. These fea- 
tures suggest that the behavior of excitons in polymer 
chains can be very close to that of an ideal semiconduc- 
tor quantum wire system. 

In this paper, we study excitons in a polymer chain 
coupled to a single electromagnetic mode in a cavity. 
The resulting mixed states of excitons and cavity photons 
can be describes in terms of polariton quasiparticles. We 
will focus on the properties of polaritons in the presence 
of a deformation potential interaction with the lattice. 
We will consider a modified Su-Schrieffer-Heeger (SSH) 5 
model, describing the propagation of excitons in the poly- 
mer—, with an additional term that takes into account the 
coupling with a single cavity mode. Excitons are mod- 
eled as excitation of two level systems localized within 
the unit cells of the polymer chain, and a variational 
mean field approach is used to calculate the properties of 
the ground state of the systems at zero temperature as a 
function of the density of polaritons. A similar approach 
was used to study the transition between a polariton Bosc 
Einstein Condensate (BEC) to a laser-like behavior for 
polaritons in a cavity<£ The coupling to the lattice by the 
deformation potential adds new features to the ground 
state properties of polaritons. We will show that a self- 
trapping^ of polaritons occurs at a threshold value for 
the polariton density. This mechanism could give rise to 
a BEC of polaritons which localize spontaneously with- 
out the need of external traps or strain fields. Com- 
pared to the semiclassical treatment we have used in the 
case without cavity^ the description in terms of cavity- 
polaritons give a clearer picture of the physics of the self- 



trapping process, and introduces new features related to 
the strong exciton-cavity photon coupling. Although our 
model contains a single cavity mode, we expect the result 
to be relevant also in the description of organic systems in 
planar cavity geometries, for which interesting interplays 
of the exciton-cavity and exciton- LO phonon dynamics 
have been predicted^ Exciton-polaritons are one of the 
most promising candidates for the realization of BECs in 
condensed matter systemsiii Exciton-polaritons lasing 12 
and matter-based parametric amplifiers^ are additional 
important applications of these quasiparticles. Exciton- 
polaritons in organic systems are particularly interesting 
due to the big excitonic oscillator strength and to their 
strong coupling to phonons, which gives rise to strong 
optical nonlinearities4^ Evidence of polaritonic effects 
in a single polydiacetylene chains has been recently re- 
ported^ We will focus on semiconductor polymer chains 
with a non-degenerate ground state. A well known exam- 
ple in this class of materials is polydiacetyleneA&. How- 
ever, we will keep our theory general in such a way that 
it can be extended to polymers with similar properties. 

The paper is organized as follows: in Sect.^we intro- 
duce the model. In Sec. IIIII we illustrate the variational 
approach used to calculate the ground state properties. 
In Sec. lIVI we compare the variational results to an exact 
calculation with two sites in order to establish the valid- 
ity of our approach. Sect. [V] provides the details of the 
numerical energy minimization procedure used to find 
the distribution of polaritons in the chain and the total 
energy of the systems. The results on the self-trapping 
phase transition as well as analytical results that can be 
obtained is some limits are described in Sec. lVIl Sec. lVIII 
analyzes the role of the excitonic spontaneous emission 
and of the finite Q-factor of the cavity. Conclusions are 
in Sec. lVnTl 



II. THE MODEL 

The model consists of the ID system of excitons cou- 
pled to lattice deformations as well as to a single cavity 
mode of the electromagnetic field. The Hamiltonian can 
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be written as 

H = Hd + Hssh • (1) 

The first term corresponds to the Dicke modelii of an 
ensemble of two level systems coupled to a single electro- 
magnetic mode 

H D = loJc + ~{Blc + c*B n ) +ujxY. B t B n ■ (2) 

n n 

where ujx is the exciton energy. We will use h = 1 
throughout the paper. B^, B n are operators of creation, 
annihilation of excitons in a singlet spin state and c\ c 
are creation, annihilation operators for the cavity pho- 
tons of energy lu c . Each site n of the model represents 
a single monomer of the polymer chain. The parameter 
g indicates the exciton-cavity dipolar coupling constant. 
In the Dicke model, the atoms do not have translational 
degrees of freedom and there is no direct transfer of opti- 
cal excitation from one atom to the other. Here, we need 
to add two additional features to the Dicke model: (i) the 
excitations can hop from site to site and move along the 
backbone of the polymer chain, and, (ii), the hopping of 
the excitation depends of the relative position of the sites 
in the lattice, which can move around their equilibrium 
position. This can be represented by an the excitation 
transfer of the SSH form as 

Hssh = +J2j( u n+i ~ u n ) 2 

n n 

- ^ t n +x,n{B n+x B n + B\B n+ i) . (3) 

n 

The SSH model was used to describe the electronic 
transport in polyacetylene chains^, and was extended to 
the case of excitonic transport in polydiacetylene^. In 
Eq. ©, M, C, u n and p n are the mass, elastic constant, 

I 



E = cj c X 2 + - \ 
Jo 



where y>(y), p(v) and u(v) are also continuous functions. 
The prime indicates the derivative with respect to the 
variable v. In deriving Eq. 10 we have assumed that 
|a(z/)| > \fi{v)\ along the chain, which corresponds to 
the condition of negative detuning between the exciton 
resonance and the cavity mode. Notice that the total 
polariton number 

N P = X 2 + J \(3\ 2 dv = X 2 + ^ - i J ^l^dv (7) 

is a conserved quantity. In order to find the ground state 
of the system at a fixed polariton density we perform a 



total displacement and momentum of the nth site of the 
chain. The hopping term is i n +i.„ = t — ~f(u n+1 — u n ), 
where 7 is related to the exciton-phonon deformation po- 
tential D = 27a. a is the site separation in the tight- 
binding chain. The value of to is determined by the ex- 
citon effective mass m as tn = „ 1 • > . 

u 2ma z 



III. ENERGY MINIMIZATION 

We consider a polariton trial wave function which is a 
product of a coherent state for photons and a BCS state 
for excitons as — 

\X,a,(3,<p) = \\)l[(a n \0) n + e i ^[3 n \l) n ) , (4) 

n 

where A, a ni /3 n , and ip n are variational parameters, and 
I A) represents a coherent state for the cavity. The co- 
efficients a n and (3 n are subject to the single-occupancy 
constraint 

\a n \ 2 + \(3 n \ 2 = 1 . (5) 
|0) is the vacuum state of the cavity mode, and |0)„ and 
|l) n denote the ground and the excited state of the two 
level system at the site n. We assume A real and we 
fix the phase tp n to make a n and (5 n real. Due to the 
hopping, the variational coefficients will depend on the 
index n in the general case. 

We can interpolate continuously the wavefunction by 
transforming the discrete sum Y2n=i m a continuous in- 
tegral J Q N dv. In this way, we can define the optical polar- 
ization Tp{v) = 2a(v)P(v), where a and f3 are continuous 
functions of v. We can express the total energy of the 
system as the expectation value of the Hamiltonian in 
Eq. Q with the trial wavefunction in Eq. @. In the 
continuum form, the total energy can be expressed as 



dv , (6) 

I 

variational minimization of (H — fiNp) with respect to 
the functions ip, tp, u and with respect to the constant 
A. From the condition 5(E — fiNp) = we obtain the 



\- + Cu' 2 + i (V' 2 + ip 2 if' 2 ) - t ijj 2 + 1U '^ 2 + g Xip cos tp- oj X \Jl-ip 2 
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system of equations 



to{<p"ilf + tp'^ ) + sin 93 = , 
Cu"+1^' = 0, 

— oX cos (p , = , 

(lu c — |Lt)A + — / ijj cos yjfiz/ = . 



(8a) 
(8b) 

(8c) 
(8d) 



From Eqs. I)8a|) and iJBJl we can see immediately that, 
in the assumption a(3 > 0, and g > we can take the 
solution ip — 7r, since this variable has no constraint. 
That makes the global phase of the optical polarization 
constant along the chain and out of phase with respect 
to the cavity field. In order to have a closed equation for 
■0 we can integrate the equation for u in Eq. I)8b|l and 
substitute in Eq. lj%c|) . By direct integration of Eq. 
we can write 



7 

2C 



h/f 



aA 



(9) 



where A is a dimensionless constant of integration. Wc 
choose A = 0, which implies that the total length of 
the polymer is not fixed. The force constant C can be 
expressed in terms of the sound velocity S as C = 
Finally, the system of equation in Eq. (JSJ can be rewritten 
in the form 

- 2i V - t Q r - X ^ -gX + { ^j=4r = ' ( 10a ) 

v 1 - V 

(lu c - fi)X - I / ipdv = , (10b) 



where the coefficient of the cubic term \ = A ® IS -i ■ This 
system of coupled equation will be solved numerically in 
Sec. E| 



IV. TWO SITES MODEL 

In order to check the validity of the BCS trial wave 
function in Eq. J3J and the accuracy of the variational 
approach, we compare in this section the exact solution 
of the problem with a result obtained with the varia- 
tional approach for two lattice sites. The exact solution 
is obtained by writing the wave function in the form 

# = j]o\00Np)+r ll \0lNp~l)+r l2 \l0Np-l)+r l3 \llN P -2) . 

where the first term in the right-hand side corresponds to 
the state with zero excitons and Np photons, the second 
and third terms describe the states with one exciton and 
Np — 1 photons, and the last term is a state with two 
excitons and Np — 2 photons. Here, we assume that the 
number of polaritons in the system Np > 1. In this 



N P 


Eq, exact 


E , trial WF 


A 

4 




3.089 


20 


16.928 


16.972 


100 


87.837 


87.857 


200 


177.012 


177.027 


1000 


893.528 


893.521 



TABLE I: Ground state energies calculated exactly and using 
the trial wave functions. The calculations were made using 



the following set of parameters: 
g = 0.2£ 9 . 



section we will also consider that the hopping ti2 as a 
fixed parameter. Using this form for the wave function in 
the Schrodinger equation with the Hamiltonian in Eq. 
we obtain a set of equations for r/f. 



irjo = lo c N p j]o + - ^/Np{rji + r] 2 ) 



ir\\ = £iVi + '^ivNpm + \/Np — lr) 3 ) - ti 2 i] 2 



im = £iV2 + ^(yNpm + V n p - l7 ?3) - turn 



ii) 3 = (lo c (Np - 2) + 2lo x ) % + ^\/Np - 1(771 + 772) , 



where e\ = uo c {Np — 1) +ujx- The energy spectrum and, 
in particular, the ground state energy are found by direct 
diagonalization. 

Using the trial wave function approach described in 
the previous section, we find that 

(a, A, 0, ip\H\a, A, /?, <p) = uj c X 2 + Wx(/3? + Pi) - 

2 

-2fi 2 ai/3ia2/3 2 cos(y>i - ip 2 ) + gX }^ ajpj cos (fj . (11) 



Without a loss of generality, the followin g substitutions 
are made in Eq. ifTTjl: ip t = n, X = N P — fi\ — f3 2 , 
a i = V 1 ~ Pf ■ The resulting expression depends only 
on Pi. The minimum value of (a, A,/3, cp\H\a, X, j3, ip) is 
considered as the ground state energy. Table [I] shows the 
comparison between the exact and the variational results. 
The ground state energy was calculated for different val- 
ues of the photon number in the cavity. The ground state 
energies calculated with the variational approach are in 
good agreement with the exact calculation. As expected, 
the exact ground state energy is slightly smaller than the 
energy calculated using the trial wave function. This dif- 
ference can be related to correlation effects not included 
in the trial wave function. Also notice that the difference 
between the ground state energy calculated with the two 
approaches decreases at a larger number of the photons. 
In fact, the correlation effects between the excitons and 
the cavity photons are expected to disappear for a pho- 
ton number much larger that the excitation number in 
the system. 
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V. NUMERICAL PROCEDURE 

Once established the validity of the variational ap- 
proach, we can study the general case of a long chain. In 
order to study on the same ground the inter-cell hopping 
and the presence of the cavity field we solve numerically 
the system of Eqs. (fill)) . We use the steepest descent 
method of functional minimization^ This method is ef- 
ficient in solving numerically Gross-Pitaevskii equations 
for BECsjSI which are similar to our nonlinear equations. 
The method consists of projecting an initial trial state 
onto the minimum of an effective energy Ti. by propa- 
gating the state in an imaginary time. We start from 
an initial field parameter A(r = 0) and a trial function 
iP(t = 0), then A(r)and V ; ( T ) are evaluated in terms of 
the equations 
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FIG. 1: Polarization ip for 5/t = 0.01, g/to = 0.01, x/to = 6, 
pp = 3 x 10~ 3 (dotted line), pp = 6 x 10" 3 (dashed line), 
pp = 6 x 10 -2 (dot-dashed line), and pp = 10 (solid line). 



dtp(T) _ sn 



and 



dr 5ip n {r) 



a\(r) _ sn 



(12) 



dr 8\{t) 



(13) 



where 8 indicates a constrained derivative that preserves 
normalization. Eqs. I|12|) and (|13fl define a trajectory 
in the parameter space for the optical polarization and 
the field parameter A. At each step we move a little 
bit down the gradient — =^ and — =p The end product 
of the iteration is the self-consistent minimization of the 
energy. The time dependence is just a label for different 
configurations. In practice we chose a step At and iterate 
the equations 



ip(r + At) w iP(t) - At 



sn 

64>(t) 



A(t + At) w A(t) - At 



sn 
sx(r) 



(14) 



(15) 



by normalizing ip and A to the total number of polaritons 
Np at each iteration. The time step At controls the rate 
of convergence. The system is described using 100 points 
and periodic boundary conditions and one parameter for 
the cavity photon field. As a test we compared our nu- 
merical calculations with some analytical limit cases. For 
the trial initial ijj we used a random values on each site 
and also a form corresponding to an analytical limit that 
will be discussed in the next section. The number of it- 
eration depends on the convergence rate and the choice 
of the initial trial function. Typically, we used 10 5 — 10 6 
iterations. 



VI. RESULTS 

A. Low polariton density 

The results of the numerical solution for the function 
ip are shown in Fig. ^ for the case x — 6£o< At very low 
excitation densities pp = ^p, where N is the number of 
sites, polaritons behave as simple bosons. In that limit 
the self-trapping effect is absent since the nonlinear at- 
tractive potential is proportional to the local polariton 
density and can be neglected in this low density limit. 
Mathematically, this limit can be described by approx- 



imating 



ip in the last term of Eq. (jTU|) . The 



cubic term ip 3 is not strong enough to give rise to the 
self-trapping effect in this limit. 

In the absence of the self-attractive term the distribu- 
tion of the optical polarization is homogeneous. In this 
low excitation limit the chemical potential fi is simply 
given by 



VP 



9'- 



(16) 



where lo' x = lux — 2to (energy at the bottom of the ex- 
citonic band), S = lo' x — lu c is the optical detuning, and 
g' = gVN- Notice that in this limit the chemical poten- 
tial does not depend on the coefficient of the cubic term 
in Eq. I|10|l . Moreover in this limit the chemical poten- 
tial corresponds to the energy of the lowest polariton at 
fc = 0. 



B. Intermediate density: self-trapping 

By increasing the number of polaritons we observe a 
critical polariton number at which a symmetry-breaking 
transition occurs. A similar symmetry- breaking transi- 
tion was found in the case of BEC with an attractive 
nonlinear interactional Our numerical calculations for 
the dependence of /i on the polariton density pp = ^f- 




FIG. 2: Polariton chemical potential as a function of the po- 
lariton density for S/to = 0.01, g/ta = 0.01, x/to = (solid 
line), x/to = 4.5 (dashed line), and x/to = 6 (dot-dashed 
line) . 



are shown in Fig. for several values of the cubic 
term x- When the exciton number increases following 
the polariton number, a self-trapping occurs, caused by 
the polariton-lattice coupling (dashed line in Fig. 
There is a drop in the chemical potential at the point 
of symmetry-breaking (Fig. |5J due to the self-trapping. 
The breaking of the spatial homogeneity gives a neg- 
ative contribution to the energy of the system. Also, 
this makes exciton-like states more energetically favor- 
able than photon-like states, which results in a sharp 
decrease of the photonic component for nonzero x, as 
shown in Fig. The self-trapping effect brings a dis- 
continuity both in the chemical potential fi, and in the 
photon (exciton) density as a function of the polariton 
density. Since the effective attractive potential depends 
on both and Xi the point of symmetry breaking for 
higher x corresponds to smaller values of the polariton 
density p P = ^ . 

Starting from the small excitation limit, we can expand 
— ^ — <x ip + hip 3 . The second term introduces an ef- 



fective repulsion related to the intrinsic fermionic nature 
of the excitons and decreases the cubic term x by " X 2 ~ M . 
Some analytical forms for the solution of Eq. (|10afl with 
fixed exciton number can be written in terms of Ellipti- 
cal functions. In particular, in the case when the photon 
number is much smaller than the exciton number, we can 
neglect the fourth term in Eq. (|lUa|) . since as seen from 
Eq. (EH, A is inversely proportional to cu c — fj,. In this 
limit the solution can be written in the form 



-0 = Asech[F(j/ — v )] , 



where 



Y = 



(X-^)N P 

4in 



.4 



YN F 



(17) 
(18) 
(19) 




FIG. 3: Photon density as a function of polariton density for 
X/to = (solid line), x/to = 4.5 (dashed line), and x/to = 6 
(dot-dashed line). 



and the chemical potential is determined by the equation 



fi = u>x - 2i 



16tn 



(20) 



C. Saturation 

By further increasing the polariton density, the inter- 
nal fermionic structure of the excitons gives rise to a 
hard-core repulsion term, which leads to the saturation 
of the exciton states, making the polarization distribu- 
tion broader and flatter at the top (dot-dashed line in 
Fig. [IJ. When the saturation spreads over the whole 
chain length, the polarization distribution (solid line in 
Fig.QJ becomes homogeneous again. We expect a discon- 
tinuity of fi again at this point as a consequence of the 
disappearance of the gradient of the polarization function 
(kinetic term) , which gives a positive contribution to the 
energy. With the further increase of the exciton density 
the hopping effect is reduced due to the blocking. Fig. 0] 

shows half of the exciton density - ^ du as a function 
of pp, which reaches 1 in the limit of large excitation 
density regardless the value of Xi corresponding to half 
filling of the exciton band. The half filling maximizes the 
polarization and hence minimizes the dipole interaction 
energy between the excitons and the cavity photons. In 
this saturation regime the polaritons become photon-like, 
since an added excitation mainly contributes to the cav- 
ity mode, and thus the chemical potential approaches lo c 
(see also Fig.0). A gradient in the density of polarization 
produces a force on the ions according to Eq. 1)8 b[) ■ The 
force is stronger at the edges of the saturation region of 
the \i/j\ 2 distribution as seen in Fig. EJ and is positive to 
the left from the center of the symmetry-breaking point 
and negative to the right. This reduces of the total length 
of the chain due to the interaction with the electromag- 
netic cavity mode. 
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FIG. 4: Half of exciton density as a function of the polariton 
density for x/to = (solid line), xAo = 4.5 (dashed line), 
and x/to = 6 (dot-dashed line). 
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FIG. 5: Force as a function of the site index for xAo 

pex — 

line) . 



6 x 10~ 3 (dashed line), p ex = 6 x 10~ 2 (dot-dashed 



VII. HOMOGENEOUS DEFORMATION AND 
ROLE OF DAMPING 



The lattice deformation in the presence of optical exci- 
tations can be understood by analyzing the total energy 
of the system in the assumption that the lattice defor- 
mation is homogeneous. As seen in the previous section, 
this assumption is justified for a systems in the satura- 
tion regime or for a polariton density below the critical 
value for self-trapping. The total energy of the lattice is 
modified by the presence of polaritons. The homogeneous 
deformation allows us to obtain an analytical expression 
for the total energy of the system, and to analyze the ef- 
fect of the finite linewidth of polaritons. We will include 
both the spontaneous emission of the excitons and the 
finite Q-factor of the cavity mode. Instead of using the 
variational approach of Sec. lIIII we will solve directly the 
equations of motion for the excitons and cavity photons. 



We start by using the Fourier transform 



ikna j 

e b k 



(21) 



where N is the number of lattice sites, k is the wave 
vector, and a is the lattice separation, to define the op- 
erator bk as the annihilation operator of an exciton with 
a wave vector k. The Hamiltonian in Eq. can then 
be rewritten as 



H = J2 



vl c 

2M + 2 {Un+1 



+ 



+ J2"kb{b k + J2f(k,k')b{b k , 

k k,k f 

+to c c)c + -r^bo + —cbl , 
where lu^ = u>x — 2£o cos(fca) and 

n 

^ — ika^ina(k-k') _|_ ^ik' a ^ina(k' — fe) 



(22) 



(23) 



We consider a finite-length chain with periodic bound- 
ary conditions. The equilibrium lattice displacement is 
homogeneous and will be indicated by u = u n+ \ — u n . 
Notice that in this limit f(k, k') is diagonal 



f{k,k') = 27 cas(ka)udk,k> 



(24) 



Eqs. H24fl and 11' I'll show that in this homogeneous case 
there is no mixing of polariton with different k vector. 
Since the cavity couples to the k — exciton mode only, 
the polariton modes at k ^ are completely decoupled 
and do not enter in the dynamics. Taking into account 
Eq. (|24p. the excitonic part in H reduces to uio(u)blbo, 
where u)q(u) = lux — 2to + 2~/u is the energy of the exciton 
at k = renormalized by the lattice deformation poten- 
tial. We also add a term that describes the pumping of 
the cavity mode by an external field and represent the 
full system by 



H = J2 



2M 2 



+uj (u)blbo + -w(c^b + b\c) + u c c^c + 



kE q 



). 



(25) 



where E represents the electric field of the external 
pump and k is the coupling between the external pump 
and the cavity mode. The operators B n are fermionic on 
the same site, but commute for different sites. Therefore, 
the commutation relation for bo reads 



N 



(26) 
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Since we will assume here that there are no excitations 
of phonon modes at k ^ 0, we have that ^ B^B n 



x 



b\bk ~ 6q6o- The equations for the expectation values 

k 

of the polarization (bo) = p, exciton density (6q&o) = n x 
and cavity photon operator (c) = A can be obtained using 
the standard factorization scheme 22 and read 



lo c X + —p 



kE q 



-iuj c t 



(27) 



p = —i 



nx 



Wo(1 _^ + ^ A (l_^ 
V N w 2 V N ' 



-Tp, (28) 



-i^(l-^)(p*A-pA*)-2rn, (29) 



where we have introduced the spontaneous emission rate 
of the excitons T and the damping of the cavity mode a c , 
due to the finite Q-factor of the cavity. The equation for 
A can be explicitly integrated and gives 



■ x 9 P 
i\ = — - — 

2 a c 



kE { 

2a~ r 



g — iwct 



(30) 



This expression can be used in Eqs. I|28|l and (|29|l which 
give for the steady state 



■ g' kE 



2 2 



£-0 ( 1 _ £H2L\ 
a c V 1 N I 



U> (1 



2nx 
N 



iT 



■ n' 2 



2n x ) 
N 



(31) 



where we have defined p = pe~ 

The u-dependent total potential energy U(u) of the 
system can be obtained by substituting the solution for 
p, A, and nx which depend on u, into the initial Hamil- 
tonian l25|) . At weak excitation we obtain 



NC 

U{u) ~ —^-u 2 +uj (u)nx(u) 



+UJ C 



2a 



2T 



nxiu) + lu c I 



(32) 



where / 



and 



n x (u) - \p\ 2 = 



{5 + 2 7 m) 2 + r 2 



(33) 



Notice that nx has a resonant behavior as a function of 
u, which is a consequence of the deformation potential 
shift in the excitonic energy. 

In the absence of light the potential energy U depends 
quadratically on u with a minimum at u = 0. When 
the laser is switched on, this dependence changes due to 
the presence of optical excitations in the system. Due 
to the resonant behavior in Eq. (|33ll . it is energetically 
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FIG. 6: Effective lattice potential as a function of displace- 
ment. The figure has been obtained for the following set of 



parameters: E g 



2t = 2.282 eV, D 



S = 2.5-KTm/s, a = 1.5 nm, M = 2.110" 21 g. 



= 27a = 6.1 eV, 
lo = 2.282 eV, 

and C — S 2 M/a 2 = 82eV/a 2 . The excitation corresponds to 
p ex ~ 0.1. etc = g' = 1 meV, and lo c = 2.242 eV. T = 400 /ueV 
and r = 700 /ieV for the upper and lower curve, respectively. 
The curves are displaced for clarity. 



favorable to have u ^ if this reduces the total energy 
of the system. Figure shows the potential energy for 
some values of parameters. The parameters used in these 
calculations are typical of polydiacetylene chains. We 
observe that the potential U has a parabolic dependence 
with a superposed Lorentian contribution due to the opti- 
cal excitation. This Lorentian contribution can be either 
positive or negative depending on the relative strength 
of the exciton-cavity coupling, spontaneous emission rate 
and cavity damping. The actual behavior can be explic- 
itly calculated using Eq. I|32l) . In the figure we consider 
two particular sets of parameters providing a minimum 
of the energy that corresponds to contraction and expan- 
sion of the lattice. 



VIII. CONCLUSION 

In conclusions, we have investigated the optically- 
induced lattice strain in a single polymer chain in a cav- 
ity. We have extended the excitonic SSH model, describ- 
ing the exciton propagation in the chain, to include the 
effect of the cavity electromagnetic field. Using a po- 
lariton picture, we have obtained a system of integro- 
diffcrcntial nonlinear equations for the spatial distribu- 
tion of the optical excitations in the chain. We find so- 
lutions describing a self-trapping of the polaritons which 
saturate when the excitation density is increased. The 
chemical potential of the polaritons shows a discontinuity 
at a threshold value of the polariton density. This criti- 
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cal density corresponds to the onset of the self-trapping. 
The self-trapping causes a sharp increase in the excitonic 
component and decrease in the photonic component of 
the polariton wavefunction. We have also considered the 
role of the finite radiative recombination rate of the ex- 
citons and the finite Q-factor of the cavity. These can 
be studied in a direct way in the case of a homogeneous 
strain field. We have found that both a contraction and 
expansion of the lattice are possible, depending of the 
relative strength of the exciton-cavity coupling, radiative 



recombination, and cavity Q-factor. 
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